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Abstract 

We present an algorithm for the asymmetric traveling salesman problem on instances which 
0^ . satisfy the triangle inequality. Like several existing algorithms, it achieves approximation ratio 

' O(logn). Unlike previous algorithms, it uses randomized rounding. 

(N 

1 Introduction 

Let y be a set of n vertices and let c : V xV ^ M-l- be a cost function. We assume the triangle in- 
equality: Cjj < Ci^k + Ck,j for all vertices i,j, k. The asymmetric traveling salesman problem (ATSP) 
^ . is to solve min^r J2vev '^v,'k(v) over all cyclic permutations vr on y. A subgraph { {v, Tr{v)) : v € V } 

Q I is called a tour; we seek a minimum cost tour. 

c/3 , We will use the following standard notation. For [/ C y, define 

5+{U) = {{v,w) : V £U, w ^U}, 
5~{U) = {{w,v) : V £U, w ^U}. 

^ ' For any vector x € R^j^^^ and F CV xV, we use the notation x{F) = YleeF^e- 

I The Held-Karp linear programming relaxation of ATSP is as follows. 

0^ ' min c^x 

s.t. x{6-{{v})) = x{6+{{v})) Vv£V 



O 



x{6~^iu)) > 1 ytD^ucv 

X > 



X 

5_j ' By standard shortcutting arguments, we may assume that < 1 for all v, w, and that 

. P^. ! x{5+{{v})) = 1 for all v. 

Several polynomial-time algorithms [3] [SJ pp. 125] [Ij [2] are known for computing a tour 
whose cost is at most a factor O(logn) larger than the optimum. In addition, several proofs [10\ [8] 
are known showing that the integrality gap of the Held-Karp relaxation is O(logn). This note 
provides another such algorithm and another such proof of the integrality gap. 



2 The Algorithm 

The algorithm proceeds in two steps. In the first step, we round the fractional solution using a 
simple randomized rounding schema to obtain nearly-balanced graph. In the second step, we solve 
the patch up problem to make the graph Eulerian. The algorithm succeeds in returning a connected 
Eulerian subgraph of small cost with high probability. 
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2.1 Constructing a nearly balanced graph 



Let X be any feasible solution to this linear program. Since x is balanced at every vertex, this 
implies that x is Eulerian, i.e., x{S^{U)) = x{5^{U)) for all U QV. So x is a fractional solution for 
which all cuts are perfectly balanced. We now use x to construct an integral solution z for which 
all cuts are nearly-balanced, i.e., for each cut U V, < 2. Moreover, the cost of z is at 

most O(logn) • c^x. 

The first observation is that x has equivalent cut values to an undirected graph. Formally, for 
U CV, define 6{U) = {{v,w} : veU,w<^U}. For y G M^^^ and FC{^), let y{F) = Y^eepVe- 

Claim 1. Since x is Eulerian, there exists y € such that y{5{U)) = x{5^{U)) for all U (^V. 

Proof. Define y{v,w} = {xv,w + Xw^v)/"^ for all v,w. Then 

~t~ Xnj^y 1 

as required. □ 



vim) = E t '^"''^ = \{x{5+m+x{5-{u))) = x{6+m, 



We now apply a random sampling result of Karger [5]. For convenience, we reprove it here 
in our notation. For any undirected graph with minimum cut value c, Karger [5j shows that the 
number of cuts of value at most ac is less than (2^) • This result applies to the graph induced by y 
and hence, by Claim [H also to x: 

|{ [/ : / [/ C x{6+{U)) <a}\ < n^". (1) 

To round x, we must first scale it so that its minimum cut value is large. Let G be the directed, 
weighted, multigraph obtained from x by taking K := 100 Inn parallel copies of each edge, each of 
the same weight as in x. Let q be the value of the i*^ cut, ordered such that K < ci < C2 < ■ ■ ■ ■ We 
will construct a directed multigraph H by taking each edge of G with probability proportional to 
its weight. The expected number of edges chosen by H in the i^^ cut is Cj. Let pi be the probability 
that the actual number of edges chosen in the i^^ cut diverges from its expectation by more than 
an e fraction. By a Chernoff bound, pi < 2e~^ '^'Z^. 

We will ensure that no cut diverges significantly from its expectation by choosing e appropri- 
ately and applying a union bound. Define e = y^l/lO. Since Ci > K = 100 Inn, we have pi < n~'^ 
for all i. For the small cuts, we use the bound 

2n2 



= 0{l/n). (2) 



i=l 



For the large cuts, we use a different bound. Eq. ([T]) implies that c„2a > aK. Letting i = n?", we 
have Ci > i^rini/(21nn), and hence pi < Thus 



^Pi < Y.'^^'^ < / ^"^^^^^ = 0(l/n). (3) 
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Combining Eq. ([2]) and Eq. ([3]) shows that with probabihty 1 — 0(l/n), no cut in H diverges from 
its expectation by more than an e fraction. The expected cost of H is K ■ c^x, so a Chernoff bound 
again implies that the cost of H is O(logn) • c^x with high probabihty. 

Let z € Z^f^^^ be the vector giving the total weight of the edges in the multigraph H. Assuming 
that no cut in H diverges significantly from its expectation, we have 

The last inequality follows because e < 1/3. Thus z is a nearly-balanced graph with high probability. 



2.2 Patching Up 

We now make z Eulerian by "patching it up" with another graph w. That is, we seek another 
vector w € l}^^^ such that z + w is connected and Eulerian — an integral, feasible solution to the 
Held-Karp relaxation of ATSP. Indeed, we show that a subgraph of z can be used to patch up z. 

Consider the transshipment problem on V where each vertex v has demand b{v) := z{6^{v)) — 
z{5~{v)). Hoffman's circulation theorem [71 Corollary 11. 2f] implies that there exists a subgraph of 
z giving a feasible, integral transshipment for these demands iff the capacity of each cut is at least 
its demand: 

z{6-m > J2b{v) = z{6+{U))-z{6-{U)). 

This inequality is implied by Eq. (j4]), so the desired transshipment w exists, and its cost is at most 
c^z. Thus z + w gives a connected, Eulerian graph of cost at most 2c^z, which is O(logn) • c^x, 
as argued above. By shortcutting, we obtain a tour of no worse cost. If x is an optimum solution 
of the linear program then the resulting tour is at most a factor O(logn) larger than the optimum 
tour. Consequently, the integrality gap of this linear program is at most O(logn). 



3 Tight Example 

We now show that the analysis of the algorithm given above is tight to within constant factors — 
we give an example where we must choose K to be r2(logn) in the first step of the algorithm. This 
condition is necessary not only to ensure the nearly-balanced condition but also to ensure that H 
is connected. 

Consider any extreme point x such that Xq < | for every arc a ^ A and let E be the support 
of x. Such extreme points exist of arbitrarily large size [T]. Using the fact that \E\ < 3n — 2 where 
n = \V\, we obtain that there in an independent set of vertices Vi of size at least ^. Since Xa < | for 
each a S yl and x{6~ (v)) + x{5'^ (v)) = 2, the probability that u is a not an isolated vertex in H is at 
most 1 — 2y7r for each v £ Vi. Since Vi is an independent set, these events are independent. Hence, 
the probability that none of the vertices in Vi is an isolated vertex in H is at most (1 — 2y7r)^- In 
order for H to be connected with constant probability, we must take K = 17 (log n). 
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